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ABSTRACT 
a computer programs GRNDMC and ARBNMC numerically determine vlf/1f mode 


conversion coefficients and mode sums for an earth-iunosphere waveguide which 

is inhomogeneous along the dtrestion of propagation. GRNDMC generates values 
of the vertical electric field, Eas at the ground as produced by a ground 
based vertical dipole source. ARBNMC calculates at an arbitrary height within 
the guide, all three electric field components (Be, Ex and Ey) generated by an 
electric dipole of arbitrary orientation and height within the guide. The 
ARBNMC program treats air-to-air, ground-to-air or air-to-ground vlf/1f 
problems. It will also handle ground-to-ground propagation but is much less 


economical to run than GRNDMC. 
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I. INTRODUCTION 


This outline is designed to help a knowledgeable user run the NOSC 
“simplified” mode conversion computer programs, GRNDMC and ARBNMC. They are 
modified and improved versions of those programs originally developed and 
described in references (1) and (2), respectively. The improvements include a 
considerable savings in computer core storage (more than a factor of three in 
the case of ARBNMC), an increase of the number of waveguide slabs from 25 to 
essentially an unlimited number, better formatting of the input data, 
especially in the reduction of from three to two cards per mode, an increase 
in the number of allowable modes from 5 to 20, and the number of modes to be 


considered may vary from one ionospheric slab to another. 


The computer programs described in references 1 and 2 were developed for 
calculating vif/1f field strengths in the earth-ionosphere waveguide when 
allowance must be made for horizontal inhomogeneity in the direction of 
propagation. They, GRNDMC and ARBNMC, are particularly relevant to the 
problem of propagating across the sunrise and sunset terminators and for 
propagation in an artificially disturbed environment. All of these programs 
are based upon a slab model, assume waveguide invariance in the direction 
normal to the great circle path between transmitter and receiver, and neglect 
reflections resulting from inhomogeneity along the direction of propagation. 
The field calculations, principally through waveguide modal constant inputs, 


allow for vertical inhomogeneity as well as anisotropy of the ionosphere. 


In the GRNDMC program, field strength calculations or mode sums are 
generated for the vertical electric field, E,» at the ground produced by a 
ground based vertical dipole. ARBNMC calculates all electric field components 
E,, E, and Ey at any receiver height within the guide (x-z is the plane of 
propagation) generated by electric dipole exciters of arbitrary orientation 
located at any height within the guide. Thus air-to-air, ground-to-air, air- 
to-ground, and ground-to-ground vl1f/1f propagation in a _ horizontally 


inhomogeneous waveguide channel may be treated using the ARBNMC program. 


Familiarity is assumed with references 4 and 5 which describe Fortran 
programs for obtaining waveguide mode constants and the excitation factors for 
electric dipoles of arbitrary orientation located at any height within the 
earth-ionosphere waveguide. Crucial outputs from these programs are the 


ground efgenangles and four independent quantities from which a tensor array 
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of nine excitation factors relating to end-on, broadside or vertical dipole 


excitation of E,, E, and Ey may be determined. These quantities for each mode 


and slab serve as input to the GRNDMC and ARBNMC programs. These quantities 
are obtained from the programs described in references 4 and 5 by setting the 
variable NPUNCH = l. 


Principal outputs of the present programs are mode conversion 
coefficients (in a generalized sense) and mode sum plots as a function of 
distance from the transmitter. Since the mode conversion coefficients are 
independent of the location of the horizontal inhomogeneity relative to the 
transmitter, provision is made for moving the inhomogeneity in increments and 
plotting mode sums for the incremented distances (this option is useful only 
if the ground conductivity and geomagnetic orientation may be regarded as 


constant over the path). 


In section II a general description of the program options is given. 
Section [II summarizes the relevant formulae. A description of the program 
input, output and operating procedures is given in section IV. Results are 
given in sections VI and VII. Appendix 1 gives a FORTRAN listing of GRNDMC 
while Appendix 2 gives a FORTRAN listing of ARBNMC. 


in litre ar ie ys ea oni entrant reeds eet yee 


i oat. 


II. GENERAL PROGRAM DESCRIPTIONS 


at aa 


| The mode conversion programs, GRNDMC and ARBNMC, provide approximate 
By methods for numerically determining mode conversion coefficients and mode sums 


for a slab model of the earth-ionosphere waveguide which is inhomogeneous 


along the direction of propagation (i.e., along the great circle path between 
| : transmitter and receiver). Invariance normal to the transmitter-receiver line 
| is assumed in the calculations and reflections due to the _ horizontal 
inhomogeneity are neglected. The computational procedure used in these 
programs is referred to as “simplified” in that a full wave fields program for 
“height gain” terms is not required. The associated height gain integrals are 
instead performed analytically as contrasted to the more accurate full wave 
integration schemes used in other less simplified propagation models, such as 


described in reference (3). 


Mode conversion programs are particularly relevant for propagation across 
the sunrise (or sunset) terminator and for propagation in an artificially 
a disturbed environment. The field strength calculations allow for vertical 


inhomogeneity as well as anisotropy of the ionosphere. 


Two distinct options are available with the present GRNDMC and ARBNMC 
programs. One option is for field calculations (amplitude and phase) as a 


function of range for a fixed location of the horizontal inhomogeneity. The i 


tae. 


second allows for fteld calculations at two distinct receiving points along i¢ 
the same great circle path as a function of position of the horizontal 
: inhomogeneity. An important limitation of this second option tis that the 


calculations are useful only if the ground conductivity and the geomagnetic 


$ parameters change little over the path. | 


F| GRNDMC (ground based propagation via mode conversion) is a modified 
¥ version of the program described in reference 1. It generates values of the 
vertical electric field, E,, at the ground as produced by a ground based 
vertical dipole source. The program output consists of mode conversion 
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coefficients and a field strength plot as a function of transmitter receiver 
distance or as a function of location of the horizontal inhomogeneity relative 


to the transmitter for fixed point transmissions. 


ARBNMC (for airborne propagation via mode conversion) is a modified 
version of the program described in reference (2). ARBNMC differs from GRNDMC 
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to the extent that it can be used to calculate all electric field components 


j 
4 E,, Ey, and Ey for any receiver height within the guide (x-z is the plane of 
a propagation) for electric dipole exciters of arbitrary orientation located at 
| any height within the guide. Thus air-to-air, ground~to-air, air-to-ground, 


t and ground-to-ground vif/lf propagation in a horizontally inhomogeneous 
waveguide channel may be treated using the ARBNMC program. The principal 

3 outputs of this program are mode conversion coefficients and mode sum plots as 
| a function of distance from the transmitter for the three electric field 
components for four orientations of the electric dipole exciter. The 

transmitter and receiver must be within the earth curvature dominated portion 

of the guide (below the ionosphere) but otherwise their altitude is 

arbitrary. Since the mode conversion coefficients are independent of the 

location of the horizontal inhomogeneity relative to the transmitter, 

provision is made for moving the inhomogeneity in increments and plotting mode 
7 sums for the incremented distances (again, this option is useful only if the 
| ground conductivity and geomagnetic orientation may be regarded as constant 


over the path). 
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III. SUMMARY OF THE EQUATIONS: 


In the propagation of vlf and if terrestrial radio waves to great 
distances, the waves are confined within the space between the earth and the 
fonosphere. This space acts as a waveguide, and the “waveguide concept” is 


applicable for characterizing the propagated fields as a function of distance. 


The waveguide mode method obtains the full wave solution for a waveguide 
that has the following characteristics: (1) arbitrary electron and ion 
density distributions and collision frequency (with height) and (2) a lower 
boundary which is a smooth homogeneous earth characterized by an adjustable 
surface conductivity and dielectric constant. This method also allows for 
earth curvature, ionospheric inhomogeneity, and anisotropy (resulting from the 


earth's magnetic field). 


The energy within the waveguide is considered to be partitioned among a 
series of modes. Each mode represents a resonant condition, t.e., for a 
discrete set of angles of incidence of the waves on the fonosphere, resonance 
occurs and energy will propagate away from the source. The complex angles (0) 
for which this occurs are called the eigenangles (or “modes"). These are 
obtained using the “full-wave" procedures described in references 4 and 5 by 


solving the determinantal equation (i.e., the modal equation): 


F (0) = [R,(0) RCO) - 1[ = 0 (1) 
Ry gq 6) By g ©) 

where Ra (9) = (2) 
Ry gq 69) yg © 


is the complex fonospheric reflection coefficient matrix looking up into the 


fonosphere from height "d" and 


RW Peat napik Se Ke te 5 eh A CTIEET ws niga: ENA AE GID BIOES O SR Ie . 2 ” Bt tee ne See 


ae 18 a ALEC PE TERE YR NES Mee! : 


iB d (9) e 
: R, (©) = (3) 
1 0 iB a (9) 


| is the complex reflection matrix looking down from height "d" towards the 
ground. 

The notation I for the R's and R's denotes vertical polarization while 

i the notation 1, denotes horizontal polarization. The first subscript on the 

R's refers to the polarization of the incident wave while the second applies 


to the polarization of the reflected wave. 


The individual terms of equations (2) and (3) are: 


wy = the ratio of the reflected field in the plane of i 
incidence to the incident field in the same plane. | 
! Ry = the ratio of the reflected field perpendicular to the ) 
* plane of incidence to the incident field perpendicular to ; 
: the plane of incidence. 
Wy = the ratio of the reflected field perpendicular to the 
’ plane of incidence to the incident field in the plane of 
incidence. 
Ry = the ratio of the reflected field in the plane of 


incidence to the incident field perpendicular to the 


plane of incidence. 


The ionospheric reflection matrix, Ry (equation 2) at height, d, may be 
obtained (from references 4 or 5) by numerical integration of differential 
equations given by reference 6. The differential equations are integrated by 
a Runge-Kutta technique starting at some height above which negligible 
7 reflection is assumed to take place. The initial condition for the 
E. integration, i.e.,the starting value of R, is taken as the value of R for a 

{ sharply-bounded ionosphere at the top of the given electron density and 
collision frequency profiles. The method for obtaining this starting solution 
| 
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is described in reference 7. The term R, may be calculated (as in references 


4 or 5) by the methods described in reference 8 in terms of solutions to 


ih aes « 
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| Stoke's equation and their derivatives. That is, the ground reflection 
coefficient matrix R, , as given in equation (3), is determined in terms of 


independent solutions, h, and hj to Stokes’ equation 


Z where the functions hy and hy are modified Hankel functions of order 1/3 
(which are linearly related to Airy functions) as defined by the Computation 


Laboratory, Cambridge, Massachusetts, reference 9. 


The propagation geometry, upon which GRNDMC and ARBNMC are based, is 
z shown in Figure 1. In this coordinate system the direction of stratification 
is the z direction and the direction of propagation is in the x - z plane. 
; The direction of z, which is the altitude coordinate measured normal to the 
curved earth's surface, is taken positive into the tonosphere. Positive x is j 
the direction of propagation and y is normal to the plane of propagation. 
Thus, the fields exhibit no y dependence but have a dependence of x of the 
form exp(-iK sin 9 x) where K is the magnitude of the free-space propagation 


vector, 0 is the angle between the direction of the propagation vector and the 


SERRE CaO 


z direction at a point in the stratified medium where the index of refraction 

is unity. All field quantities are assumed to have an exp (iwt) dependence 

where w is the angular frequency. The radiating dipole source for the 

i => 

7 propagated fields is denoted in Figure 1 by M. The dipole is oriented within 

3 the earth-ifonosphere waveguide by the inclination angle y and azimuthal angle 

, : 

" For a vertical dipole, y = 0° and 4 = 0°. For a horizontal dipole y = 
90° and since > is the angle between the direction of the horizontal dipole 
and the direction of propagation, = 0° represents end-fire and $ = 90° 

a represents broadside launching. 

The procedure used in the mode conversion model segments the earth- 
fonosphere waveguide into M cascaded uniform waveguides (Figure 2) where the 
propagation characteristics do not change within any segment. The slabs are 
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numbered from left to right with the total number of slabs being denoted as 
M. The slab labeled NTR is the slab containing the transmitter. The modal 
equation, equation 1, is solved by the method of either references 4 or 5 for 
as many modes (eigenangles, 8,5) as desired in each waveguide segment (or 
slab). The model allows for an arbitrary number of modes on each side of any 
given boundary. At any transition region (i.e., slab boundary) energy 


scatters from a given mode into several other modes in the adjoining slab. 


Subject to the assumptions of invariance in the y-direction and the 
neglection of reflections from horizontal inhomogeneities, mode conversion 
coefficients may be computed by enforcing the boundary condition (at each slab 
interface) that all tangential field components must be continuous. Also, the 
mode conversion algorithm formulation assumes that a unit amplitude wave 
exists for each mode k in the transmitter region (slab NTR of figure 2). The 
quantities aie relate the energy from an incident mode, k, in the transmitter 
slab to j-th mode in slab p. 

The “accumulative” mode conversion coefficient aie for the p-th slab 
associated with conversion from the k-th to the j-th mode expressed in terms 
of the coefficients for the previous (p-1)-th slab may be written in the form 


(reference 10). 


J 
P PsP. 7Ps»p-l = 
: ay ry rk >» P= NTR + 1 (5) 


J 
“ Pm) (_yy.P7h a P»p~l 
A aik [ ike (x, x14 Ty NIR+ 1S <M 


where i = (-1) , K is the free-space wave number, S; is the sine of the qth 


eigenangle for the pth slab, and J is the total number of modes assumed to be 


{mportant in the total field determinations. 


A simplified explanation of equation (5) is demonstrated in figure 3. In 


the figure some important variables are: 
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a The “normalized conversion coefficient” for the p-th slab, I . 

| i> jr -1 i> 

E | 1 

tie + 7P a i+ yPoP e y yP>P- 6 

i where all terms are matrices. 7 

The “total or accumulative conversion coefficient" for the p-th slab, aP. 
: k> i> 

gi P e 7P 

z jy (2%) jt (75,) (7) 

R 

i 2 i = 

Fi -~4" oP = 

| exp[ it S; (x, x1)! 0 0 
c | “te 20 op-l 

i ° it 0 exp{-i-> Si (x > pl 0 

a Zs 2a op-l 

- 0 0 exp[-iy St i x1)! 
e°| 

= > 

: p-l 
} iy (0%, ) 
! where again all terms are matrices. The indices are j = 1, 2, 3 ...J for 
( the modes in slab p; i = 1, 2, 3 ...I for the modes in slab P-1l; and k = 1, 2, 

mm 3 ..-K to the modes in the transmitter slab NTR. 

: If tt is assumed that the transmitter is located in slab 1 (see figure 


3), the required expansions for the first 4 slabs, supposing 3 modes, are: 


[cei ee ae 
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The symbol under the various variables 


where \ is the free space wavelength. 
(i.e., a and I) denotes a matrix. 


Critical for the solution of the system of equation (5) is the evaluation 


of the integral 


ymP = f ay - GP dz (9) 


where the t denotes the adjoint and GP is a four-element column matrix of 
height gains for the y and z components of the electric and magnetic field 


strength of the kt) node in the pth slab. 


The column vector GP (z) is given by: 
Pp 

e yn (2) 

: eh?) 

Ge(z) =| yp (10) 

k hy, (2) 
Pp 

hoy l?? 


The term ay is a four-element column matrix of height gains for an appropriate 
adjoint waveguide (see reference 11). 

In the mode conversion model the height gain elements of A are assumed to 
be simply rearrangements of the elements of G and are expressed below a height 
b in the guide in terms of the modified Hankel functions of order 1/3. The 


crucial integrals in (5) are then given in reference 10 by 


j 


? Pp 2 Pp e TOPHT m e P m e P 
Te 7 68) +S) is (e eP +3 Hyp) dz (1) 


yj yk yj 


In arriving at the final form of Equation (11), the height gains below z = 0 
and above z = TOPHT have been discarded. Also, introduced is the relation: 


-az/2 
H(z) e a (12) 
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It can then be shown that in the earth curvature dominated portion of the 


guide, Hy is linearly expressible in terms of solutions to Stokes’ equation. 


This fact is also true for the e's so that if the e's and H's are assumed to 


obey Stokes' equation in the range 0 < z < TOPHT, Equation (11) becomes 


m m p 
j,k RK (s-sP) yk dq yj dq yk dq yj dq “o 
if m# p and/or j} #k 
2s" dew du” 
ee | yi,2 yj, 2 52 m \2,4b 
gg + GN + accer + A], (13) 


if m = p and j = k 


where ey and Hy are given in terms of modified Hankel functions of order 


1/3. Explicitly, these relationships are 


F,h,(q,) + F,h,(q,) 
ae - Le - Z a z B (4) 
y 1) Io 2%) Io 


Fh, (q,) + F,h,(q,) 
Jo. 2 L - : = 2 § (15) 
y 3%) I a) Io 


where “b" = TOPHT in equations (13), (14) and (15). The terms Fy, Fo, F3 and 


F, are given by equation (16) through (19). 
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n 1/3 
) ak 2 _ 2 yi/2 
F, = - {H,(a,) - 4 = Ap ay - hy (a,)} 
2 
n 1/3 
Fy = Ca) - £27 2 ada, 
s 
t 1/3 
F,-- {hy @)-1 Gp a,” - 8°y"/? nea} 
1/3 
. ak 2 2.1/2 
-2/3 
a,7 qm  (? -2 (r2)) 
2/3 
HC) = nya) +5 A (% "a hy(q) 5 j= 1,2 


2 2 
rn, 1- gq fn z) 
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where: 


C = cosine of the angle of incidence at height h 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


h = the height at which the modified refractive index is taken to unity 


K = the free space wave number 

e/e, = dielectric constant of the ground 
o = the ground conductivity 

w = the circular radio frequency 


a = the earth's radius 


The subscript on q represents the value of z at which q is evaluated. 


example, q, means that (20) is to be evaluated for z = 0 and q, means that 


2 


(20) is to be evaluated at z * b. Similarly, the subscript on n° represents 


the value of z for which (22) is to be evaluated. 
17 
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For 


Again, the functions h, and h, are modified Hankel functions of order 1/3 
as defined by the computation Laboratory at Cambridge, Massachusetts, and the 


primes on these quantities denote derivatives with respect to the argument. 


The variable "b" in equations (20) and (22) is defined as the effective 
ionospheric height. The value to be assigned to this variable has been 
determined (reference 1) to be that height level for which the conductivity 


parameter, w is: 


r’? 


Wy, = 2.5 * 10° sec7! (24) 


This value of "b"” is the reference height (h'), as defined by Wait (reference 
12) in terms of exponential profiles. In the present programs, the value of 


“b" is allowed to vary from slab to slab. 


The column vector G,.(z) of equation (10) is uniquely defined by 
normalizing the y component of the rf magnetic field to unity at the ground 
and by introducing the proper amount of TE wave at the ground which from mode 


theory is given in reference 10 by: 


Ps (1 + ,Ry) CL ~ Ry Ry) E (1+ RY) Ry 48, ss 
h = = a 
y semi ad adh bs J eae Ges ane ee ae ee 
d=0 d=0 


This value of "f" is to be substituted into equation (15). 


Inputs to the GRNDMC and the ARBNMC programs for each slab "P" and each 
mode "“j" in the slab are the ground eigenangles (0';)'s and the T,'s as 
defined below. The T,'s are readily obtainable from either of the mode 
finding programs of references 4 or 5, as are the ground eigenangles. These 


“T" quantities are: 


POOREST CBS 


1/2 ae = 
sh? + B)2a - OR, GR) 
Hs aL an aes 
F'CO 5) Ry Py d 
1/2 — 2 = 
gh/24a + R720 - WR, GR) 
1° Ce an 
F'(,) iB Dp, F 
: sta + RG. + RD AR es 
3 F(,) Dip : 
R 
ae |} 
oe) eo (29) 
4 Wt ¢ 


where S is the sine of the eigenangle at the reference height “h" and F'(O5) 
{s the derivative of the mode equation, equation 1, evaluated at the 
eigenangle, 0; The R and R's represent, respectively, elements of the 
reflection matrix looking into the fonosphere and towards the ground from the 
level z = d. Consistent with the usual notation, the first subscript refers 
to the polarization of the incident wave and the second subscript refers to 


the polarization of the reflected wave and 


2 
Dio 7 CFR (ay) + Fohy(ay) (Fah, (a4) + Fyho(a4)) (31) 
Doo = (Fshy (ag) + Fyhglag))” (32) 


again F,, Fo, F3, and F, are given by equations (16) through (19). 
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The mode conversion programs, as they are now programmed, require that 
the height variable “d" be at the ground so that in equations (26) through 
(32), 22d 0. 


The function f of equation (25), which is the ratio ey to hy at the 
ground is computed from the T's of equations (26) through (29) as: 


f-— or +--+ (33) 
1 3 4 


t, = DL, T » tT * Doo Ty » 137 Dio T° (34) 


In terms of the above quantities the excitation tensor elements are given as: 


Field Component > Exciter? 
Vertical 
A= QD) re ren 
i — Gane 


The columns relate to excitation of the electric field components E,, E, and 
Ey and the rows apply to excitation by a vertical dipole Q,), a horizontal 
dipole end-on Og) and a horizontal dipole broadside (Ag)- Recall the 
geometry of the situation is such that z is taken positive into the 
fonosphere, that positive x is the direction of propagation and that y is 


normal to the plane of propagation. 

Both the modal excitation factor and the modal height gain functions are 
needed in computing electric field strengths. The excitation factors of 
equation (35) must be supplemented with the height gains which are defined as 


following: 
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£,(2) = exp (24) (F,h,(q) + Fyh,(a))/(Fh,(4,) + Fohy(a,)) (36) 


ae, 
f,(2) "kOe 


£,(2) = (Fyh,(q) + Fzhy(q))£/(Fgh,(4,) + Fyho(a,)) 


where f, is the height gain for the vertical electric field E,, fy is the 
height gain for the horizontal electric field component E, and f3 is for the 


electric field component Ey which is normal to the plane of propagation. 


In terms of the excitation factors, height gains and generalized mode 
conversion coefficients, the electric field components EP in the p-th slab may 


be written, as a function of receiver range (Pp), as follows: 


NTR NTR NTR NTR NTR 


E = i (A f (z,.) cos Y +4 f. (z,,) sin Y cos $¢ 
5 
n [si ( fa] 2 Ink 1k T 2nk 2k T 


NTR .NTR NTR - 
+2 f (27) sin Y sin 6) fk (zp) e 


ik (NTR 
3nk *3k "i 


k - L)p ; 


p = NTR 


+ DR, IRR Rey 8 CRORES mg 


gP = —__2___ i) Gat Ga cos Y + at, Go sin Y cos $ 


lnk 1k 2nk 2k 


n  [sin(x/a)]!/? 5 k 


P 
S 
NTR NTR 2 NTR, p..p 7 p j 
+ Asn fan (27) sin Y sin 6) (6), + (1 - 63,) 8° /8))E, (Zp) je NTR 
k 
NTR Pp 
~ixtS +s, (P -x_) - 9) 
wg ik "k *yTr+1 7 $4 p aten (40) 


The receiver altitude is Zp and the transmitter altitude zp. The final k 
subscript on the A's and f's denotes mode indices whereas the index n takes on 
the values 1, 2 and 3. Consistent with the previous definition n = 1+E,» 
n = 2>E, and n = 3+E,. Also, the notation “NTR” refers to the transmitter 


slab. The constant Q is 
Q = 0.03248k/Y F (41) 


with the free space wavenumber, k, in inverse km and F the frequency in kHz. 


The symbol Sj represents the Kronecker delta. That is 
Ss = 1 i= j (42) 
fig 948 

The angles y and 4 determine the orientation of the electric dipole source 


| 
relative to the x, y, zZ2 coordinate system as shown in Figure l. When 
executing the GRNDMC program, only the vertical E. field is computed, the 
transmitter and receiver heights are both zero and Y = $ =0. In both GRNDMC | 
and ARBNMC field strength amplitude is expressed in dB above a microvolt per 
meter for a one kilowatt radiator and phase in degrees relative to free 


space. 
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il IV. DESCRIPTION OF PROGRAM EXECUTION 


| A. General Comments 


To handle horizontal inhomogeneities, the ionosphere is divided into a 
j series of vertical slabs, as described in section III. These slabs are 
§ 3 labeled 1, 2, - . . , M as shown in Figure 2 and the boundaries between the 
| slabs have coordinates Xj, %3, + + + , Xye For each slab and for each mode 
the ground eigenangle (9') and the T,'s defined by Eqs. (26) through (29) mst 

be provided. Note that the transmitter slab is denoted as “NTR”. 


Several variables which are functions of the earth's magnetic field or 
are related to the ground conditions must first be identified. The 
geomagnetic field is specified by three variables. "AZIM" is the clockwise 
angle between magnetic north and the horizontal propagation direction in 
degrees east of north. "“CODIP” is the complement of the magnetic dip angle in 
degrees. The magnetic equator is specified by "CODIP = 90". "MAGFLD" is the 
magnetic field intensity in webers per square meter. The ground conditions 
are specified by two variables. "SIGMA" is the conductivity in mhos per meter 
and "EPSR" is the relative dielectric constant. The mode finding programs, 
references 4 or 5, punch the values of these variables on the first output 


card of each "slab”. This card is denoted as the RFACMSET card. The 


vartables on this card are identified as: 


R: The slab reference distance in megameters. 


F: The propagation frequency, in kHz, “FREQ”. 


% 
¥ 
9 : The geomagnetic azimuth in degrees, "AZIM". 
Ss C: The geomagnetic codip in degrees, “CODIP”. 
‘al 
Pt M: The magnitude of the geomagnetic field in weber/square meter, 
' “MAGFLD". 
| S: The earth's conductivity in mhos/meter, “SIGMA”. 
3 f E: The relative dielectric constant of the earth, "“EPSR". 
; T: The height of the top of the slab, “TOPHT"”. t 
| 
E | 
{ 
| 
{ 
q 23 
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B. Description of Input 


All input to the mode conversion program is given in a data deck on the 
standard input unit. Listings of sample input showing data deck setup are 
given in example I and example II. Example I illustrates the input cards for 
the option IPLTOP = 1, while example II shows the input for the option 
IPLTOP = 2. 


There are two parts to the input. The first part is read in by means of 
a Fortran NAMELIST input format. The first card of each set of input must 
contain a blank in column 1 and &DATUM in columns 2-7. This is followed by at 
least one blank and then data items separated by commas. The data items have 
the following forms: (all cards must have a blank in column 1) 

variable name = constant, 

or 

array name = set of constants, (all separated by commas), 

The second part of the input follows the NAMELIST input. The first input 
for this part is an identification card. It has u> to 40 columns of alphanu- 
meric information and follows the control card, "DATA". The information on 
the identification card is used to label output plots. Following the 
identification card a series of punched cards (obtained from the programs 
described in references 4 or 5 with NPUNCH = 1) is input for each slab of the 
modeled earth-ionosphere waveguide. The first card is denoted as the RFACMSET 
and gives the value for R, FREQ, AZIM, CODIP, MAGFLD, SIGMA, EPSR and TOPHT. 
Next there are two cards per mode. The first of these contains the complex 


eigenangle at the ground (i.e., O', and 0';) and values for the variables T, 


r 
(Real), T, (Imaginary), Ty (Real) and Ty (Imaginary). See equations (26) and 
(27). The second card contains the eigenangles at the ground (Ory and 0',) 
(duplicate input) and Tz (Real), T3 (Imaginary), T, (Real) and Ty 
(Imaginary). See equations (28) and (29). If the number of modes in any 
given slab is “NRMODE”, there will be 2*NRMODE cards for each slab. The 
2*NRMODE cards for slab No. 2 follow those for slab No. 1 and so on up to slab 
number M. Although ordering of modes is not critical, they are ordered 
according to their real parts. (The mode with the largest real part is called 


mode 1.) 


The following variables and arrays may be specified in the NAMELIST 
input: 
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IPLTOP ~ Plotting option flag. If IPLTOP = 1, two plots (field amplitude 
in dB above a uv/m for 1 kw radiated power versus transmitter ~terminator 
distance for two receiver positions) are produced for each set of input. If 
IPLTOP = 2, a plot (field amplitude in dB above a uv/m for 1 kw versus 


distance from transmitter) is produced. 


XVAL ~- The numerical values (in km) of the X's (see Figure 2) which are 
the horizontal positions of the boundaries between adjacent slabs. Note that 
XVAL can be negative and that it is dimensioned for 50. Figure 2 illustrates 
that the XVAL distances are taken relative to the transmitter. The slab 
number is identified by ISLAB and a particular slab is identified by 
XVAL(ISLAB). That is, XVAL(3) = X3- Note that the terminator line is taken 
to be the boundary XVAL(2). 


RCDOPT ~- Option to use the RFACMSET-card values (i.e., R) for XVAL(C ISLAB) 
values. If RCDOPT is set equal to one and IPLTOP is set equal to two, then 
the values used for XVAL(ISLAB) are taken from the RFACMSET card rather than 


from NAMELIST. 


RHOMAX - Maximum horizontal distance in km at which field strengths are 


desired. 


RHOMIN - Minimum horizontal distance in km at which field strengths are 


desired. 


DELRHO - Horizontal increment in km for which successive field strengths 


are computed. 


NTMAX - Number of times the transmitter - terminator separation is 


incremented. (NTMAX = 1 for IPLTOP = 2.) 


DELTAX - Distance in km by which transmitter - terminator separation is 
incremented. (DELTAX = 0.0 for IPLTOP = 2.) 


GAMMA or INCL = Dipole orientation angle relative to z axis (see 
Figure 1). Note that GAMMA is dimensioned for 4. 
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relative to X axis (see 


PHI or THETA - Dipole orientation angle 
Figure 1). Note that PHI is dimensioned for 4. 
NRP - Number of GAMMA-PHI pairs up to 4. 


TALT - Transmitter altitude in km. 


RALT - Receiver altitude in km. 


ICOMP ~ The electric field components to be computed. ICOMP = 1 is E,, 
ICOMP = 2 is E, and ICOMP = 3 is Ey: 
ICMPMX - The maximum number of components which will be computed. If 


equals 1 only EZ. If equals 2 only E, and E,. If equals 3, E,, E, and Ey 


will be computed. Note that if only E, is wanted, ICMPMX still must be set 


equal to 2. If only Ey is required, ICMPMX still must be set equal to 3. 


TOPHT -~ Height above ground in km above which height gains are 
discarded. Note that a TOPHT for each slab is required and dimensioned for 
50. If the value of TOPHT(ISLAB) is zero then the values of TOPHT used in the 
internal calculations within the program are taken from the RFACMSET cards. 


H - Is the height in km at which the modified refractive index is unity 
(n2 = 1 - 2/a(H - z)). a is the earth's radius. H is also the height to 


which the eigenangles are referred. 


TH - Option which, when set to one converts the SIN(9') to SIN(0},). The 
values of 0' are the eigenangles at the ground while the values QO, are the 


eigenangles referenced at the height H. The equation is 
sin(9,) = (1 - H/a) sin(0'). 
This relation is important in computations of the excitation factors. 


INIFLG - Printing option flag. INTFLG must be set to 1 if printout of 
height gain integrals and the "normalized" conversion coefficients are 


required (i.e., Teer, yPoPm} and iy of equation 6). 
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4 IPRNTA - Printing option flag. IPRNTA mist be set to 1 if printout of 
the “accumulated” or total mode conversion coefficients is required 


i (i.e., aP of equation 7). 
1 


NPRINT ~— Printing option flag. NPRINT equal to 1 gives a brief summary 


Pinata acai hi Mad lee ER 


i of slab input data including slab number and the number of modes for each 
| slab. Printed out are the values of all the variables contained on the 
RFACMSET cards. NPRINT = 3 gives the same printed output as for NPRINT = 1, 
except that the values of the ground eigenangle and the values of the T's of 


equations (26) through (29) are also given for each slab. 


ITXRX - Printing option flag. ITXRX must be set to 1 for debugging 


purposes. 


NAPLOT - Plotting option flag. NAPLOT must be set to 1 if plots of 


amplitude vs. distance are required. 


NPPLOT - Plotting option flag. NPPLOT must be set to 1 if plots of phase 


3 vs. distance are required. 
SIZEX ~ Length of x-axis in inches. 


SIZEY ~ Length of y-axis in inches (for amplitude). 


ot ge eee ere 


‘ XMIN - The value (in km) of the x-coordinate at the beginning of the 


x~axis. 


YMIN - The value (in dB) of the y-coordinate at the beginning of the 


y-axis. 


XINC - The change in the x-coordinate value (in km) between successive 


tee 


a) labeled tic marks. 


YINC - The change in the y-coordinate value (in dB) between successive 


labeled tic marks. 
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XTIC - The distance (in inches) between tic marks along the x-axis. 


YTIC - The distance (in inches) between tic marks along the y-axis (for 


amplitude). 
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NTICX - The repeat cycle for placing coordinate values at tic marks on 


the x-axis. 


NTICY - The repeat cycle for placing coordinate values at tic marks on 


the y-axis (for amplitude). 
SIZEP - Length of y-axis in inches (for phase). 


PHSMIN - The value (in degrees) of the y~-coordinate at the beginning of 


the y-axis. 


PHSINC - The change in the y-coordinate value (in degrees) between 


successive labeled tic marks. 


sie commit a ee ee 


PTIC - The distance (in inches) between tic marks along the y-axis (for 


phase). 


NTICP ~ The repeat cycle for placing coordinate values at tic marks on i 


the y-axis (for phase). 


FACT - The ratio of the desired plot size to the normal plot size. For 


Se RRR Nyse ee 


example, if FACT = 2.0, all subsequent pen movements will be twice their 
normal size. When FACT is reset to 1.0, all plotting is normal size. 


Note: For NTICX, NTICY and NTICP: 


= 1 - causes values to be placed at every tic mark 


| = 2 - causes values to be placed at every second tic mark, etc. 


z = 0 - suppresses all coordinate values 


The end of the NAMELIST input is signaled by &END. 


Initial values of the NAMELIST variables are presented in Table 1. 
28 
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TABLE 1 


NAMELIST VARIABLES AND INITIAL VALUES 


NAME NAME UNITS 


*ICOMP 
* I CMPMX 
LTXRX 
INTFLG 
IPRNTA 
IPLTOP 
NRP 


SIZEX inches 
SIZEY inches 
XMIN km 
YMIN db 
XINC km 
YINC db 
XTIC inches 
NTMAX YTIC inches 
NPRINT NTICX 

RHOMAX 0.0 NTICY 

RHOMIN 0.0 RCDOPT 

DELRHO 0.0 IH 

DELTAX 0.0 H 

*GAMMA 4*0.0 degrees *RALT 

*PHI 4*0 .0 degrees *TALT 

* INCL 4*0.0 degrees SIZEP inches 
*THETA 4*0.0 degrees PHSMIN degrees 
XVAL 50*0.0 kn PHSINC degrees 
TOPHT 50*0 .0 PTIC inches 
NAPLOT NTICP -- 
NPPLOT 

FACT 


~ —- ~~ O00 0 ww & 


Note: Variables denoted with asterisk (*) apply only to the ARBNMC program and 
not to the GRNDMC program. 
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3 C. Special Rules 


Some special rules mst be applied to the selection of some of the 
NAMELIST variables. 


(1) Consider the following relationship as related to plotting: 


| 4X - increment of x-scale in km/inch 
{ 


; AX 


= XINC { km } 
Ei (NTICX) * (XTIC) 


inch 


v] 
to 
2 
8 
x 
a 


Example (a). Given: XINC 
NTICX 
XTIC = 1 inch 


w 
i) 


2000 km 


oe 1s tnoh 1000.0 km/inch 


4X 


Also: 


AY - increment of y-scale in db/inch 


YINC db 
: on (NTICY)(YTIC) { inch } 
| 
5 
: Example (b). Given: YINC = 30 dB 
‘ NTICY = 2 
3 YTIC = 0.5 inch 
4 


30 dB 
AY 2 (0.5) inch 30 dB/inch 
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Further: 


Example (c). Given: XINC = 1000 km 
NTICX = 1 
XTIC = 1 inch 


1000 km 


ax =< (i) inch 


= 1000 km/inch 


Note (1): 


Example (c) shows that, if values of the axis are to be labeled at every 
tic mark, then NTICX = 1. Also, if XTIC is 1 inch, then AX = XINC. 
Therefore, in this case XINC indicates the value of the increment of the 
x-scale in km/inch. In the same manner, YINC would indicate the value of the 


increment of the y-scale in dB/inch. 


(2) Consider the following relationships: 


XMAX - XMIN 
XMAX - XMIN 
oes Dettax 7? 


Note (2): 


XMAX is not a NAMELIST variable but it must be considered when deciding 
on the values for the other NAMELIST variables. 


Example (d). Given: DELTAX = 200.0 km 
AX 1000.0 km/inch 
XMIN = -2000 km 
“XMAX" = 5000 km 
SIZEX = ? inches 
NTMAX = ? inches 
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XMAX - XMIN 
Ax 


SIZEX = 


5000. ~(-2000) 


100070 =z 7 inches 


SIZEX = 


XMAX - XMIN 
NIMAX = ——DeLTax * } 


5000. -(-2000) 


200 + 1 = 36 


NTMAX * 
(3) The following rules must be used when picking values of certain NAMELIST 
variables: 
For “IPLTOP" = 1, Set: 
XVAL(1) = -9999.9 


XVAL(2) = XMIN 


For "IPLTOP" = 2, Set: 


XVAL(1) = 0.0 
XMIN = 0.0 
NTMAX = 1 

DELTAX = 0.0 


RHOMIN = DELRHO 
RHOMAX < "KMAX" 


(4) XVAL distances are taken relative to the transmitter (NTR). The 
terminator line is XVAL(2). 


(5) The maximum value possible for SIZEY or SIZEP is 8 inches. 


(6) For the option IPLTOP equal to 1, the criteria for modeling “sunset” or 
“sunrise” at the receiver are shown in figure 4a through d. Figure 4a, for 
sunset illustrates that the daytime mode constants (smaller value of TOPHT) 


are assumed for slab No. 1, while the nighttime mode constants are assumed for 
32 


ah sean age ck ee 


ry 


ee er eee 


ie der as Con 


le ee LEED 
; c 


~ & 


slab No. 4. The mode conversion programs allow for incrementing the XVAL 
distances (which are distances relative to the transmitter) by the distance 
DELTAX. Figure 4b shows the resulting mode sum value in dB as a function of 
the distance between the terminator line (XVAL(2)) and the transmitter. 
Figures 4c and 4d illustrate the input and resulting mode sum curves for the 
“sunrise” condition. In this case the mode constants for nighttime are 
assumed for slab No. 1 while the daytime mode constants are assumed for slab 
4. 


(7) The mode conversion programs, as they are now formulated, require that 
the height “d" of equations 1, 2, 3 and 25 through 32 be set to the ground 
level (i.e., d = 0). This requires that when the T's of equations 26 through 
32 are obtained from the program described in reference 4, the variable "D” in 
that program must be set to zero. Also, if the T's are obtained from the 
program described in reference 5, the variable "GRDFLG" in that program must 


be set to l. 


(8) The maximum number of data points (i.e., distances) for which fields can 


be computed is 402. 


(9) The maximum number of allowable modes is 20. 


(10) The maximum number of allowable slabs is 50. This value may be 
increased by increasing the dimensions of the variables TOPHT and XVAL. 


D. Program Layout 


The subroutines used in the GRNDMC and the ARBNMC mode conversion 


computer programs are shown in Table II. 


® 
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TABLE II 


Program Subroutines 


GRNDMC 


MAIN 
CLINEQ 
AXISM 
MAGANG 
MDHNKL 
MCSTEP 
HTIMTL 
*MCFLD (for GRNDMC) 
MCFLD2 
*MCPLTS (for GRNDMC) 
MCPLT2 


ARBNMC 


MAIN 
CLINEQ 

AXISM 

MAGANG 

MDHNKL 

MCSTEP 

HTIMTL 

*MCFLD (for ARBNMC) 
ACCUMA 

*MCPLTS (for ARBNMC) 
HTGAIN 

CURVE 


The flow of the programs GRNDMC and ARBNMC is shown in Figure 5. Note 
that GRNDMC utilizes two temporary disc files and that ARBNMC utilizes four 


temporary disc files. 


The subroutines denoted with the asterisk indicate that these subroutines 


are not identical in the two programs, 
34 
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| SUBROUTINE HTGAIN (Z, H, HFLAG) 


Z is dimensioned for 2. 2Z(1) is set equal to the transmitter height TALT 
| and 2(2) is set equal to the receiver height RALT. The height gain functions 
fi» fo and f3 defined by Eqs. (36), (37) and (38) respectively are computed 
for the transmitter and receiver heights. H is the height where the modified 
| index of refraction is one. IH = 0 indicates H = 0 and IH = 1 indicates H # 0 


I Tee 


in the subroutine calculations. 


SUBROUTINE HTINTL (IFLG, P, INTFLG) 


HTINTL calculates the height integrals defined by Eq. (6). NORM (1P»P) 
- is an array of 20 by 20 which contains all combinations of modal integrals for 
= the slab p. Also CAPI (rP>P-!) 45 an array of 20 by 20 which contains all 
combinations of modal height gain integrals for the slab p and the previous 
: slab p-l. IFLG is a control flag set to zero in MAIN if the slab p equals 
a one. It is set to 1 if p is not equal to one. INTFLG is a printing option 
flag. It must be set to 1 if printout of NORM, CAPI and INORM (I ”) is 
desired. Note from equation 6 that the array INORM is the normalized 


conversion coefficient. 


SUBROUTINE MDHNKL (Z, Hl, H2, HIPRME, H2PRME) 
MDHNKL calculates for argument Z two independent solutions (Hl and H2) 
and their derivatives (HLPRME and H2PRME) of Stokes’ equation by methods 


described in reference 9. 


x : SUBROUTINE MAGANG (ARG, MAG, ANGLE) 


ql MAGANG converts complex number ARG to polar form with ANGLE in degrees. 


SUBROUTINE CLINEQ (A, B, X, N, NDIM, IFLAG, ERR) 


CLINEQ computes the solution of simultaneous linear equations with 


; 4 
<A] 

| 

5 complex coefficients. That is, it solves the matrix equation 
A*X=B 
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for the vector X of length N, given the matrix A of size N by N and the vector 
The A is destroyed by CLINEQ, 
IFLAG 


B of length N by Crout's L-U decomposition. 
NDIM is an integer variable which must be greater than or equal to N. 
Setting IFLAG = 1 bypasses the 
ERR is 


is an integer variable normally set to zero. 
L-U decomposition of A when solutions are required for different B's. 


a real variable computed by CLINEQ which indicates the relative errors in the 


computed solution vector X. 


SUBROUTINE MCSTEP (M) 


MCSTEP calls for CLINEQ and provides as its output the mode conversion 


coefficients for the slabs NTR, NIR + 1. . ., M for all values of NTR 


The mode conversion coefficients defined by 
“A = TOTAL CONVERSION COEFFICIENTS". 


consistent with the input data. 


equation 7 are printed out under the 


SUBROUTINE AXISM 


AXISM draws and labels the axes for the outout plots. 


SUBROUTINE CURVE (for ARBNMC) 


CURVE plots the data points on the output plots. 


SUBROUTINE MCFLD (for GRNDMC) 

MCELD called from MAIN if IPLTOP = 1 computes the vertical electric field 
at the ground produced by a ground based vertical electric dipole at RHOMIN 
and RHOMAX for transmitter-terminator distances ranging between XVAL (2) and 
NTMAX*DELTAX + XVAL (2) at intervals of DELTAX, using equations (39) and (40) 
with y = 0°, > = 0°, Z, = 0, 2, = 0 and f, = f, = f, = 1. 


SUBROUTINE MCFLD2 (for GRNDMC) 

MCFLD2 called from MAIN if IPLTOP = 2 computes the vertical electric 
field at the ground produced by a ground based vertical electric dipole from 
RHOMIN to RHOMAX at DELRHO intervals using equations (39) and (40) for a fixed 
horizontal inhomogeneity. Again, y = 0°, ¢ = 0°, Zr = Zp = 0 and f, = fo = fy 
= 1, 
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SUBROUTINE MCPLTS (ITYPE) (for GRNDMC) 


MCPLTS generates two plots (field amplitude in dB above a wv/m for 1 kw 
radiated power (or phase in degrees) versus transmitter-terminator distance 
for two receiver positions) each time it is called from MCFLD. ITYPE = 1 
gives amplitude plots. ITYPE = 2 gives phase plots. 


SUBROUTINE MCPLT2 (ITYPE) (for GRNDMC) 


MCPLT2 generates one plot (field amplitude in dB above a uv/m for 1 kw 
radiated power (or phase in degrees) versus distance from transmitter) each 
time it is called from MCFLD2. ITYPE = 1 gives amplitude plots. ITYPE = 2 
gives phase plots. 


SUBROUTINE ACCUMA (for ARBNMC) 


ACCUMA is the chief control routine of ARBNMC. The routine sets up the 
accumulated values of the "total" or accumulated conversion coefficient as a 


function of slab position. 


SUBROUTINE MCFLD (for ARBNMC) 


MCFLD when called from MAIN with IPLTOP = 1 computes the field 
components BP defined by Eqs. (39) and (40) for transmitter height Zy (TALT) 
and receiver height Zp, (RALT) for as many as four (GAMMA, PHI) pairs. GAMMA 
and PHT describe the orfentation of the electric dipole source. Calculations 
are made for ranges RHOMIN and RHOMAX for distances between the transmitter 
and the start of the horizontal inhomogeneity ranging between XVAL (2) and 
NTMAX*DELTAX + XVAL (2) at intervals of DELTAX using Eqs. (39) and (40). 
Field amplitude outputs are in dB above a uv/m for 1 kw radiated power and 


phase angles are in degrees relative to free space phase. 


MCFLD when called from MAIN with IPLTOP = 2 computes the field 
components EP defined by Eqs. (39) and (40) for transmitter height Zep (TALT) 
and receiver height Zp (RALT) for as many as four (GAMMA, PHI) pairs. 
Calculations are made for transmitter-receiver distances ranging from RHOMIN 
to RHOMAX at DELRHO intervals using Eqs. (39) and (40) for a fixed horizontal 
inhomogeneity. Field amplitude outputs are dB above a uv/m for 1 kw radiated 
power and phase angles are in degrees relative to free space phase. 
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SUBROUTINE MCPLTS (ITYPE) (for ARBNMC) 


E | MCPLTS, for IPLTOP = 1, generates six plots (three field component 
amplitudes in dB above a pv/m for 1 kw radiated power (or phase in degrees) 
| versus distance between transmitter and start of the horizontal inhomogeneity 

for two receive: ranges). As many as four (GAMMA, PHI) pairs are possible so 


that each plot can contain as many as four curves. 


q MCPLTS, for IPLTOP = 2, generates three plots (three field component 
| amplitudes in dB above a uv/m for 1 kw radiated power (or phase in degrees) 
versus transmitter receiver distance for a single location of the horizontal 
inhomogeneity). As many as four (GAMMA, PHI) pairs are possible so that each 

plot can contain as many as four curves. 


ITYPE = 1 gives amplitude plots. ITYPE = 2 gives phase plots. 


E. Description of Output 


s Sample output from program ARBNMC is shown in example III. This output 
is obtained from the input shown in example I. The output listing begins with 
a an abbreviated listing of the NAMELIST input variables. Note that the IPLTOP 
: = 2 option is requested. This is followed by printout of the control card 

“DATA” and the fdentification card, "PRESTON U.K. TO SONDRESTROM". Next comes 
i the printout of the NPRINT = 3 option giving slab number, the RFACMSET card 
values and then a listing of ©' angles and T's for all slabs and modes. The 
©' angles are the eigenangles at the ground while the T's are the complex 


quantities given by equations (26) through (29). 


z The principal output of the mode conversion program then follows. Since 
‘ INTFLG was set equal to one, values of certain parameters are printed out for 


each slab, P. These quantities are INORM (J,K) (i-e., {P 


5k of equation 6), 


which are the normalized conversion coefficients of mode K in slab P-l to mode 
% J in slab P; NORM (I,J) (i.e., uy of equation 13), which are values of the 
| integral between height gains of the I and J modes in slab P; and CAPI (I,K) 

(i.e, aie of equation 13), which are values of the integrals between 
height gains for mode K in slab P-l and mode I in slab P. Also, since IPRNTA 
was set equal to one the “total (or accumulative)” mode conversion 
coef ficifents avn of equation 7 are printed out. The tabulation represents the 
conversion from mode K in the transmitter slab to mode J in slab P. 
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Finally, the mode sum values for the vertical electric field (E,) at 
transmitter height of 0.0 km and receiver height of 6.0 km are listed in dB 
above a microvolt per meter for a1 kilowatt radiator as well as the phase in 
degrees referenced to free space propagation. In this example mode sums are 
shown for receiver locations at a distance (p) with p varying from 25 km to 
4000 km at 25 km intervals. A plot of the computed fields is shown in figure 


6 as compared with propagation data described in reference 13. 


Example IV illustrates the output obtained from the input of example 
II. In this case the IPLTOP = 1 option was requested. The first part of the 
printout shows the NAMELIST values. Because NPRINT was set to one only a 
brief summary of the input data (i.e., the RFACMSEL values) is printed after 
the "DATA" and identification cards. Also, since neither INTFLG or [PRNTA was 
set to one, no printout of the height gain integrals or conversion 
coefficients are printed out. Since IPLTOP was set equal to 1 the sample 
output shows mode sums for the three electric field components E, (ICOMP = 1), 
ey (ICOMP = 2) and Ey (ICOMP = 3) as a function of transmitter-receiver 
distance ranging from RHOMIN to RHOMAX at DELRHO intervals. It should be 
pointed out that the NAMELIST variable, ICOMP, must be entered in three 
separate calls to NAMELIST. The mode sums are listed in dB/yv/m for a one 
kilowatt radiator and the phases in degrees relative to free space. Because 
NRP = 4 in the input, there are four GAMMA-PHI pairs (i.e., four antenna 
orientations) for which the mode sums are computed. Note that TALT = 10 km 
and RALT = 10 km. Shown in figure 7 (a, b, c, d, e, f) are plots generated by 


the mode conversion program for this case. 


39 


| PI OM Te ad Sr ee ee ee F : | ee a ne ——nere A 


ME ky 1 Sty Oe Comme EO TEE ET PETS Keon! . : : vet ree a 


pS a ha as Eg oe ‘ oe Fn en CENTRE 


V. REFERENCES 


1. Pappert, R. A., and Shockey, L. R. [1975] Effective Lonospheric Height For 


| a_ Simplified Mode Conversion Model at VLF. Naval Electronics 
Laboratory Center Interim Report 761. 


4 | : 2. Pappert, R.- A. and Shockey, L. R. [1976] Simplified VLF/LF Mode Conversion 
; eo Program with Allowance for Elevated, Arbitrarily Oriented Electric 


Dipole Antennas. Naval Electronics Laboratory Center Interim Report 
77T, NTIS ADAO33412, Springfield, VA. 


3. Pappert, R- A. and Shockey, L. R. [1972} Mode Conversion Program for an 
Inhomogeneous Anisotropic Ionosphere. Naval Electronics Laboratory 


Center Interim Report 722, NTIS, AD743948, Springfield, VA. 


. 4. Sheddy, C. H., Pappert, R. A-, Gough, Y. and Moler, W. F. [1968] A Fortran 


Program for Mode Constants In An Earth~Ionosphere Waveguide. Naval 


Electronics Laboratory Center Interim Report 683. 


a 5. Morfitt, D- G. and Shellman, C. H. [1976] MODESRCH, an Improved Computer 
' Program for Obtaining ELF/VLF/LF Mode Constants. Naval Electronics 


Laboratory Center Interim Report 77T, NTIS, ADAQ32573, Springfield, 
VA. 


6. Budden, K. G., [1955] The Numerical Solution of Differential Equations 


¥ Governing Reflexion of Long Radio Waves from the Ionosphere, 

" Proceedings Royal Society (London) A227, pp. 516-537. 

2% 7. Sheddy, C. H. [1968] “A General Analytic Solution for Reflection from a 
4 Sharply Bounded Anisotropic Ionosphere,” Radio Science, V. 3, No. 8, 


3. Pappert, R. A., E- E. Gossard, and I. J. Rothmuller, {1967], "A Numerical 
Investigation of Classical Approximations Used in VLF Propagation,” 
Radio Science, V. 2, pp. 387-400. 


40 


ae 


’ 
ry 
4 


Pe Sy OEE TA THR SO ES Ke” 


Ne ee nen tananaeuataemiand 


9. Staff of the Computation Laboratory at Cambridge, Massachusetts, [1945] 
Tables of the Modified Hankel Functions as Order One-Third and of 
their Derivatives. (Harvard Univ. Press, Cambridge, MA). 


rc 


10. Pappert, R. A. and Shockey, L. R. [1974] A Simplified Mode Conversion 


Program for VLF Propagation In the Earth-Ionosphere Waveguide. Naval 
Electronics Laboratory Center Interim Report 751. 


1l. Pappert, R. A. and Smith, R. R. [1972] Orthogonality of VLF Height Gains in 
the Earth Ionosphere Waveguide. Radio Science, V. 7, pp. 275-278. 


12. Wait, J. R. and Spies, K. P. []965] "Influence of Finite Ground 
Conductivity on the Propagation of VLF Radio Waves,” Radio Science, 
NBS Journ. of Research, 69D, pp. 1359-1373. 


13. Burgess, B. R. [1972], Royal Aircraft Establishment, Farnborough, 
England, letter to Commanding Officer, Naval Ocean Systems Center, 
San Diego, CA. 


i 
% 
oy 
: 
2 
? 


VI. FIGURES 


| 
! 
i 
H 
| 
{ 
{ 


z) 


Xx 


Figure 1. Dipole M orientation within the waveguide where y 
is the inclination and ¢ the azimuthal orientation. 
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Figure 2. Mode conversion model. 
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Figure 3. Identification of the terms: “accumulative (or total)” conversion coefficient, aP 


and “normalized” conversion coefficient, TP. 
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Figure 4. SUNRISE-SUNSET geometry. 
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